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ABSTRACT 

High-redshift gamma-ray bursts have several advantages for the study of the distant universe, providing unique information about the 
structure and properties of the galaxies in which they exploded. Spectroscopic identification with large ground-based telescopes has 
improved our knowledge of the class of such distant events. We present the multi-wavelength analysis of the high-z Swift gamma-ray 
burst GRB 140515A (z = 6.327). The best estimate of the neutral hydrogen fraction of the intergalactic medium (IGM) towards the 
burst is xu] < 0.002. The spectral absorption lines detected for this event are the weakest lines ever observed in gamma-ray burst 
afterglows, suggesting that GRB 140515A exploded in a very low density environment. Its circum-burst medium is characterised by 
an average extinction (Ay ~ 0.1) that seems to be typical of z > 6 events. The observed multi-band light curves are explained either 
with a very flat injected spectrum (p = 1.7) or with a multi-component emission (p = 2.1). In the second case a long-lasting central 
engine activity is needed in order to explain the late time X-ray emission. The possible origin of GRB 140515A from a Pop III (or 
from a Pop II stars with local environment enriched by Pop III) massive star is unlikely. 

Key words. Gamma-ray burst: general - Gamma-ray burst: individual (GRB 140515A) - Galaxies: high-redshift - intergalactic 
medium 


1. Introduction 

A better understanding of the chemical enrichment and evolu¬ 
tion of the high-redshift universe is one of the fundamental goals 
of modern astrophysics. High redshift surveys have been per- 
formed by means of wide held surveys of bright quasars (e.g. 
iFan 20m or deep held analyses to iden tify distant galaxies by 
their drop-out (e.g. iBouwens et al. 2014I) . The identiheation of 
high-redshift Gamma-Ray Bursts (GRBs) add a different and 
prohtable view of the distant universe (see ISalvaterra 201.51 for 
a recent review). With respect to other probes, GRBs have many 
advantages: (i) they are detected at higher redshifts; (ii) they are 
independent on the galaxy brightness; (iii) they do not suffer of 


usual biases affecting optical/NIR surveys; (iv) they reside in av¬ 
erage cosmic regions. High-z GRBs can provide fundamental, 
and in some cases unique, information about the early stages of 
structure formation and the properties of the galaxies in which 
they blow up. For example, GRBs can be used to trace the cos¬ 
mic star formation rate (Kistler et al. 2009; Ishida et al. 2011; 
Roberson & Ellis 2012), to pinpoint high-z galaxies and explore 
their metal and dust content (Tanvir et al. 2012; Salvaterra et al. 
2013; Elliott et al. 2015), and to shed light on the re-ionization 
history (Gallerani et al. 2008; McQuinn et al. 2008), to con¬ 
strain the dark matter particle mass (de Souza et al. 2013) and 
the amount of non-Gaussianity present in the primordial den- 


1 









A. Melandri et al. 2015: GRB 140515A 


sity field (Maio et al. 2013), and to measure the level of the lo¬ 
cal inter-galactic radiation field (Inoue et al. 2010). Additionally, 
they could also provide direct and/or indirect evidences for the 
existence of the hrst, massive, metal-free stars, the so-called 
Population III stars (Campisi et al. 2011; Toma et al. 2011; Wang 
etal. 2012; Maetal. 2015). 

Since the launch of the Swift satellite (iGehrels et. al 20041) 
8 events have been identihed at redshift greater than ~ 6, 
and for 5 of them spectroscopic redshift was secured, includ¬ 
ing in the list GRB 140515A that we are discussing in this 
paper. Remarkably, some of them showed fairly bright early- 
time afterglows, e ven detectable by small robotic telesc opes 
(e.g. GRB 050904; iTagliaferri et al. 20051: l^r et al. 2006h . but 
in general the observational features of high-z events do not 
seem to differ signihcantly from those of their closer siblings, 
as clearly pointed out studying, e.g., the event with the high- 
est spectroscopically confirmed redshift, z ~ 8.2: GRB 090423 
(ISalvaterra et al. 20091; iTanvir et al. 2009h . 

In this paper we describe our observations of high-z 
GRB 140515A in Sect.l^and the result of the analyses in Sect.|3] 
A discussion is reported in Sect. |4] and the main conclusions of 
our work are summarised in Sect.|5] 

Throughout the paper, distances are computed assum¬ 
ing a ACDM-univers e with Hn = 71 kms~* Mpc~‘, - 
0.27, and Qa = 0.73 (iLarson et al. 201 ll iKomatsu et al. 20 fill . 
Magnitudes are in the AB system and errors are at Icr conhdence 
level. Raw and reduced data not explicitly reported in tables are 
available from the authors upon request. 


2. Observations 


On 2014 May 15 at 09:12:36 UT (= To), the 5wi7f/BAT triggere d 
and located the long GRB 140515A (ID’ Avanzo et al. 20141) . 
SwiftjXKY promptly detected the afterglow emission whereas 
5w//f/UVOT did not identify any credible bright optical 
candidate. A faint optical afterglow was later identihed 
by ground -based observa tions with 8m Gemini-North 
telescope dFong et. al 20 1^, the 2.5m NOT telescope 
(Ide Ug^t e Postigo et al. 20 14al) . the 2.2m GROND tele- 
scope (Graham et al. 20141) . and the 3.6m TNG telescope 
dMelandri et al. 2014a). 

Spectroscopi c observations perform ed with the Gemini- 
North telescope (C hornock et al. 2014al) and the GTC telescope 
dde Ugarte Postigo et al. 2014bl) dete^d sharp decrement in 
Hux below 8900 A, caus ed by Lycr absorption at redshift z=6.327 
dChornock et al. 2014bl) . Spectral analysis will be described in 
detail in Section 3.3.1. 


3. Results 

3.1. BAT temporal and spectral analysis 

The Swift/BAJ' data were processed with the standard Swift 
analysis software included in the NASA’s HEASARC soft¬ 
ware (HEASOFT, ver. 6.16) and the relevant latest calibration 
hies. For each GRB, we extracted mask-weighted, background- 
subtracted light curves and spectra with the batmaskwtevt and 
batbinevt tasks in FTOOLS. The mask weighted light curve 
shows a double-peaked structure. The hrst pulse started at Tq-22 
s and peaked at Tq - 18 s. It was followed by a second brighter 
pulse between To - 10 s and To -H 4 s (see Fig.[TJ. The total du¬ 
ration of the burst event in the 15 - 150 keV energy band in the 
observer frame is 790 = (23.4 + 2.1) s (90% conhdence level), 
corresponding to ~ 3.2 s in the rest-frame. 



Fig. 1. BAT mask-weighted light curve showing the count rate 
in the 15 - 150 keV energy range. 


A ht to a simple power law of the time-averaged spectrum 
from To - 22 s to To + 4 s gives a photon index E = 1.86 ± 
0.14 (x^ - 61.21, d.o.f.= 56). A power law with an exponential 
cutoff gives a moderately better ht (x^ - 54.57, d.o.f.= 55; F- 
test probability P= 98.8%). For this model the photon index is 
r = 0.99/Qgg, Fpk = 51.8/^22 0 huence in the 

15 - 150 keV band is Fbat = (6.531°'4’) x 10^’ erg/cml The 
1-sec peak Hux measured from Tq + 1.50 s in the 15 - 150 keV 
band is /pk,BAT = 0.86 ± 0.10 ph cm“^ s~'). 

We also tested the presence of a blackbody component in 
the prompt emission spectrum. We added a blackbody compo¬ 
nent to the non thermal (power-law) spectrum and we ht this 
model to the data, obtaining a photon index F = 1 and a 

blackbody temperature IcT - 1 2 . 4 / 4 ^ keV, that in the source rest 
frame corresponds to kTii = 90.3/;29 0 keV. This model provides 
an adequate ht (x^ - 54.29, d.o.f.= 54) but not a signihcant 
improvement compared to the cutoff power-law model. 

We searched for spectral evolution between the two main 
emission episodes of the prompt emission. The spectrum of the 
hrst peak (from Tq-22 s to Tq- 14 s) can be modelled as a power- 
law spectrum, with a photon index T = 2.01 /q 32- The spectrum’s 
second peak (from Tq - 14 s to Tq + 4 s) is better represented 
by a power law with an exponential cutoff (F-test probability 
P= 99.6%), with photon index T = 0.82/°®, £pk = 52.7/g° 
keV. 

The total bolometric (rest frame 1-10^ keV) isotropic 
energy, assuming the exponential cutoff model, is Fiso = 
(5.8 ± 0.6) X 10® erg at z = 6.327 and Fpk.u = 379.7/®['^ 
keV, co nsistent with the rf - E\^„ correlation within i ts 1 cr 
scatter (lAmati 20061; lAmati et al. 20081; iNava et al. 20121) . The 
isotropic luminosity is Liso = (3.6 + 0.8) x 10^^ erg s“', con- 
sistent with the Fnkrf - correlatio n within its Icr scatter 
(lYonetoku et al. 20^lNava et al. 20l3) . 

3.2. XRT temporal and spectral analysis 

The XRT began observing ~60 s after the BAT trigger, with the 
hrst 9 s in Windowed Timing (WT) mode and the remainder in 
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Fig. 2. Observed multi-wavelengths light curve of 
GRB 140515A. We report the best fits for the optical (blue) and 
X-rays (grey) band. The orange line shows the possible decay 
of the late time X-rays afterglow assuming a possible jet-break. 
Inset: The time interval of the X-ray bump in linear scale. Red 
curves represent three independent flaring episodes. 


Photon Counting (PC) mo de. We collected t he XRT data from 
the online Burst Analyzer (lEvans et al. 201 Oh and converted the 
observed [0.3-10] keV count rate into flux at 3 keV. The X-ray 
light-curve at that energy (Fig. |2]l is well described by an initial 
steep decay (ai ~ 2.9) followed by a broad bump after ~ 10^ s 
that lasted for almost 1 day. From Fig. [T] it appears clear that 
the “real” burst began at Tq^grb = To,real ~ To - 22 s. If we 
consider To,grb as the beginning of the burst, shifting backwards 
the temporal axis of Fig. |2l the initial decay becomes less steep 
(«! ~ 2.4), but it makes no difference for late times breaks and 
decay indices (see Tab[TJ. 

At later times a further steepening in the light-curve is de¬ 
tected, also confirmed by Chandra deep upper limit. If we con¬ 
sider only XRT data, the late time decay index is aiate = 3.9+0.6, 
while if we take also into account the dee p late time upper limi t 
this value appears to be at least > 2.6 (iMargutti et al. 201^ . 
Although it cannot be confirmed by the optical/NIR data, the 
late time X-ray decay index might suggest a possible jet-break 
origin. 

We performed a time-integrated spectral analysis of the X- 
ray emission from f- To = 86.4 s to f- To = 10^ s using XSPFC 
ver. 12.8.2. The best fit to the data is an absorbed power law 
model: the Galactic absorption is kept fixed to the value A! 


Gal 


2.54 X 10^° cm ^ (IWillingale et al. 20131) . the photon index i 


IS 


r = (1.79 + 0.11) and the local absorption at z = 6.327 is Nh < 
5.6 X 10^^ cm“^ (c-stat=341.1, d.o.f.= 360). Selecting photons 
from f-To = 86.4 s to f-To = 2152 s, when the spectral changes 
are minor, we obtain an upper limit on the column density of 
Nh < 8.9 X 10^2 cm-2. 

Also for the XRT spectrum we tested the presence of a pos¬ 
sible blackbody component. The resulting fit adding a thermal 
component to the power-law spectrum did not improve the ab¬ 
sorbed power law model described above. 

The total isotropic energy in the 0.3 - 30 keV rest frame en¬ 
ergy band is Fx.iso = (7.19±0.43)xl0^^erg,thus GRB 140515A 
is also consist ent with the Fx.isn - Fnk rf - Fisn correlation w ithin 
its 2 cr scatter dBernardini et al. 201:^ Margutti et al. 201^ . 


Table 1. Observed X-ray light-curve fitting results (x^jd.o.f. - 
59.31/51 = 1.16 assuming that Tq.grb = To, and ;\fVd.o.f. = 
55.89/51 = 1.09 if instead To,grb = To,real)- 


Parameter 

To.grb - To 

Tq.grb - To.reai 

tti 

2.80 ± 0.22 

2.34 + 0.15 

Tbreak 

732 + 25 s 

729 + 18 s 

0-2 

-3.87 ± 0.74 

-4.03 + 0.74 

Tpeak 

3028 ± 192 s 

2930 + 133 s 


1.03 + 0.05 

1.03 + 0.04 

Tjet-break 

> 10^ s 

> 10^ s 

folate 

3.9 + 0.6 

3.9 + 0.6 


3.3. Optical/NIR temporal analysis 

The Swift/IJNOT began observing the field of GRB 140515A 
3.7 ks after the trigger (ID’ Avanzo et al. 20141) . The afterglow 
was not detected in any of 7 UVOT Alters. This is consistent 
with the redshift reported by Chornock et al. (2014a). In or¬ 
der to provide deep upper limits, we co-added the exposures 
within the first sequence of observations (00599037000). We de¬ 
termined the count rate using a 5 arcsec circular source region 
centred at the optical afterglow position reported by Fong et al. 
(2014), and a circular background region of radius 20 arcsec po¬ 
sitioned on a blank area of sky situated near to the source posi¬ 
tion. The photometry was extracted using the UVOT tool uvot- 
source. The count rates were co nverted to magnitude s using the 
UVOT photometric zero points (iBreeveld et al. 201 ih . We used 
Heasoft software version 6.15.1 and UVOT calibration version 
20130118. 

In Tab.|2]we summarise ultraviolet and optical 3cr upper lim¬ 
its and other optical/near infrared detections of the optical after¬ 
glow. We note that at the redshift of GRB 140515A the SDSS-z 
Alter is slightly affected by the absorption of the intergalactic 
medium (IGM) since the optical depth for the Fy-a at z > 6 rises 
dramatically (i.e. Fan et al. 2006). Extrapolating the IGM prop¬ 
erties from low redshift to z = 6.327 we estimated the expected 
correction for the z' filter to be > 0.30/“;“^ mag (see Japelj et 
al. 2012 for details of the method). Despite the fact that this is 
formally a lower limit of the correction - not taking into account 
the rise in optical depth at z > 6 - we took it into account before 
converting all the observed magnitudes reported in Tab. |2] into 
flux densities. 

Although the optical light curve is sparsely sampled (Fig. EJ 
it is possible to estimate the decay index of the optical afterglow 
in the z'-band at late times. From 6 ks after the burst event the 
optical afterglow follows a power-law decay with a,' = 0.89 + 
0 . 02 . 


3.4. Optical/NIR spectral analysis 
3.4.1. GTC spectrum 


We obtained spectroscopy of t he afterglow of GRB 140515A 
with OSIRIS ( CeBa_eLaL^OOO|) at the 10 .4m Gran Telescopio 
Canarias (Ide IJgarte Postigo et al. 2014bb . The observations 
were obtained between 22:37:31 UT and 00:09:46 UT (mean 
epoch 14.184 hr after the GRB onset) with 0.6" seeing and con¬ 
sisted of 3 X 1800 s exposures. We used the R2500IVPH grism, 
which covers the range between 7330 and 10000 A at a resolu¬ 
tion of ~1600 using a 1" slit. 

The data were reduced in a standard way (bias subtraction, 
pixel-to-pixel response correction, cosmic ray removal, wave¬ 
length calibration, ID extraction, flux calibration, and com- 
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Table 2. Optical observations. Magnitudes are in AB systen 
and have not been corrected for Galactic absorption along th i 
line of sight (E(b-v) = 0.02 mag, ISchlaflv & Finkbeiner 201 11) 
References for data taken from the GCNs are: 1) Fong et a] 
2014; 2) Graham et al. 2014. 


Tmid 

Exposure 

Filter 

Mag 

Ref. 

[sec] 

[sec] 




3743.0 

841 

white 

>22.14 

UVOT 

6795.0 

393 

uvw2 

> 20.73 

UVOT 

9203.0 

1082 

uvml 

> 20.90 

UVOT 

6885.0 

1141 

uvw\ 

>21.28 

UVOT 

6982.0 

1437 

u 

>21.35 

UVOT 

3875.0 

549 

b 

> 21.13 

UVOT 

4220.0 

568 

V 

> 20.27 

UVOT 

6408.0 

480 

z' 

20.27 ±0.11 

1 

42877.1 

1500 

z! 

22.18 ±0.19 

NOT 

46462.8 

60 

z' 

22.21 ± 0.35 

X-shooter 

53650.4 

1500 

z! 

22.32 ± 0.19 

NOT 

59705.5 

1500 

zl 

22.35 ± 0.20 

NOT 

61200.0 

3000 

z' 

22.1 ± 0.1 

2 

56160.0 

1800 

J 

20.63 ±0.15 

TNG 

61200.0 

2400 

J 

20.9 ± 0.2 

2 

52344.0 

3600 

H 

20.61 ± 0.10 

TNG 

61200.0 

2400 

H 

20.9 ± 0.2 

2 



Absorption redshift [Ly a) 

Fig. 3. Fya effective optical depth in the line of sight of 
GRB 140515A compai'ed with previous GRB and QSO works. 
The coloured area shows the optical depth found by Songaila 
(2004) while grey points are measurements from Fan et al. 
(2006) with sample of quasars. 


bination of spectra) using self-made routines based on IRAF 
dTodv 19931) . The resulting combined GTC spectrum shows a 
strong continuum above ~ 8900 A, where the signal-to-noise ra¬ 
tio is ~ 20 per pixel, or ~ 40 per resolution element. 

3.4.2. X-shooter spectrum 

We observed the field of GRB 140515A with the X-shooter spec¬ 
trograph mounted at the FSO/VFT using the nodding mode with 
1x2 binning. The spectrum was acquired on 2014 May 16, start¬ 
ing at 00:42:43 UT (~ 15.5 hr after the GRB onset) and consisted 
of 2x4x600 s exposures, for a total integration time of 4800 s on 
source, covering the range between ~3000 and ~ 24000 A. The 
mid expose time is 16.3 hr (~ 0.68 d) after the GRB trigger. 
The final reduced spectrum has a signal-to-noise ratio of ~ 3 per 
pixeQ, with a seeing of ~ 0.9" (measured from combined 2D 
spectrum in the VIS and NIR arms). The flux calibration of the 
X-shooter, which is problematic in general (Kriihler et al. 2015; 
Japelj et al. 2015), is uncertain due to unavailable standard spec- 
trophotometric star in the night when the observations were done 
and because the photometric observations, which could be used 
to check the quality of calibration, have rather high errors at this 
epoch. Thus, as the photometric observations at this epoch have 
rather high errors, it is not possible to use them to reliably rescale 
the spectrum. 

3.4.3. Lyo- forest constraints on the IGM 

We analysed the ionisation state of the IGM using the Gunn & 
Peterson (1965) optical depth, defined as = - ln(7“), where 
T the average transmission in a redshift bin. Following Songaila 
& Cowie (2002) and Songaila (2004) we normalised the GTC 
spectrum (as its signal-to-noise ratio, SNR, is better than the X- 
shooter one, see section 3.4.4) by fitting a power law to the con- 

* The quoted difference in S/N between the GTC and Xshooter spec¬ 
tra is due partly to the different pixel size of the two instruments and 
partly to the better observing conditions of the GTC observation. 


tinuum, and divided it into redshift bins of 0.1 between z = 5.2 
and z = 6.3. The results are presented in Fig.[3]and in Table[3] 
We only see sky line residuals up to z ~ 5.5, above 
which we can just give detection limits based on the noise 
spectrum. Our limits are less restrictive than the ones pre¬ 
sented by Chornock et al. (2014b) due to the lower SNR, 
but show the same behaviour (Fig. [3]l. Results coming from 
both GRB 140515A and GRB 130606A (IChornocket al. 20131: 
ICastro-Tirado et al. 20151: [Hartoog et al. 201 ^ are consi stent 
with quasar measurements (ISongaila 20()4 : Fan et al. 2006h . 


3.4.4. Lycr red damping wing fitting 

We tried to fit the strongest feature seen in the spectrum (at ~ 
8900 A) to an absorption Fyman-a feature with a Voigt profile. 
Following Chornock et al. (2014b), we first computed a Voigt 
model using the same constraints, obtaining inconsistent results. 
This could be due to the fact that they do not seem to consider 
the instrumental profile, whose effect on the Ly-a feature is not 
negligible at this resolution when log(VHi) < 19. Booking at 
Fig.m we can observe the residuals of a sky line subtraction few 
angstroms blue-wards the wing, precisely at the zone crucial to 
fit a Voigt model. After a careful inspection on the 2D images of 
both GTC and X-shooter instruments, we concluded that there is 
no flux at this zone. Consequently, the wing profile is too sharp to 
get a satisfactory fit, suggesting that the absorption is dominated 
by the IGM and that the host absorption is masked. 

We then built up IGM models following the prescription of 
Miralda-Fscude (1998), fixing the lower redshift value to z = 6.0 
because the contribution to the wing shape below this redshift is 
negligible (it starts to be important closer to the host). Our best 
fit, with z = 6.3298 + 0.0004 and a fraction of neutral hydro¬ 
gen xhi < 0.002, is shown in Fig. |4] We caution that due to the 
sharpness of the wing, the few points we have because of GTC 
resolution, and the sky line next to the absorption, any formal 
constraints on these quantities would be unreliable, so the values 
should be interpreted as the most plausible estimations that we 
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Fig. 4. Left: Best IGM damping wing fit to the spectrum of GRB 
140515A. Right: Redshift evolution of the hydrogen neutral frac¬ 
tion. The dotted line shows the Gnedin & Kaurov (2014) model 
and points (see legend) the observational measurements of this 
quantity. Points with arrows are lower/upper limits. 


can obtain from the data. Moreover, especially by the fact that z 
cannot be determined by metal lines, hybrid models cannot offer 
a more accurate fit than the one showed in Fig. |4] so no con¬ 
straints on the host HI abu ndance can be derived from this event 
(for further discussion, see lMiralda-Escude 1998h . However, due 
to the sharpness of the red damping wing, it is obvious that the 
neutral hydrogen present in the IGM cannot mask neither the 
presence of a DLA nor a subDLA, as their damping wings would 
be easily identihed. Consequently, we can establish a conserva¬ 
tive upper limit of log(AHi) ^ 18.5 for the HI abundance in the 
host galaxy of GRB 140515A. As shown in Fig.|4] the fraction of 
neutral hydrogen derived from this analysis is in good agreement 
with the model by Gnedin & Kaurov (2014), and it provides a 
very relevant observational constraint. 

Last, we estimated the 3cr upper limits on the observer- 
frame equivalent width (EW) for the Si II /il260, O I /il302, 
and C II A1334. We hnd a value of 0.67 A, 1.06 A, and 1.30 
A, respectively. These estimates are a factor ~ 2 more strin¬ 
gent of what reported by Chornock et al. (2014b), resulting to 
upper limits on the gas-phase abundances of [Si/H] < -1.4, 
[O/H] < -1.1, and [C/H] < -1.0. Eurthermore, these lines are 
weaker t han the average rest-frame EW s observed for a typi¬ 
cal GRB (Ide Ugarte Postigo et al. 2012|) . In fact, the strength of 
those lines compared to the average GRB spectrum that can be 
estimated with the use of the line strength parameter (ESP, as de- 
hned in de Ugarte Postigo et al. 2012), is ESP < -3.15, < -3.89, 
and < -2.88, respectively. This means that these lines are very 
weak and that GRB 140515A exploded in a relatively low den¬ 
sity environment. However, our limits on the metals abundances 
do not allow us to put a stringent limit on the metallicity of the 
progenitor. 

3.4.5. Spectral energy distribution 

We constructed a broadband spectral energy distribution (SED) 
using the flux-calibrated X-shooter spectrum and Swift X- 
ray data (Eig. |5]l. X-shooter data are treated as outlined in 


Table 3. IGM absorption towards GRB 140515A. 


z 

T 

lim(r) 

^GP 

lim(T^" ) 

5.25 

- 

0.0594 

- 

2.82 

5.35 

0.1174 

0.0709 

2.14 

2.65 

5.45 

0.1038 

0.0767 

2.27 

2.57 

5.55 

- 

0.0739 


2.61 

5.65 

- 

0.0604 

- 

2.81 

5.75 

- 

0.0527 

- 

2.94 

5.85 

- 

0.0775 

- 

2.56 

5.95 

_ 

0.0784 

- 

2.55 

6.05 

- 

0.0614 


2.79 

6.15 

- 

0.0700 

- 

2.66 

6.25 

- 

0.0965 

- 

2.34 


Japelj et al. (2015). The spectrum was corrected for Galactic 
extinction (E(b-v) = 0.02 mag) by using the Cardelli etl 
al. (1989) extinction curve and Galactic extinction maps 
(ISchlaflv & Einkbeiner 201 11) . Regions of telluric absorption 
were masked out. Spectrum was re-binned in bins of approx¬ 
imately 50 A to increase S/N and to guarantee a comparable 
weight of the optical and X-ray SED part. The absolute flux 
calibration was hne-tuned with simultaneously obtained near- 
infrared photometric observations. The X-ray part of the SED 
was built from time integrated observations obtained between 5- 
60 ks and its mean epoch was interpolated to the mean epoch of 
the X-shooter observations. 

The SED fitting w as carried ou t with the spectral htting 
package XSPECvl2.8 (lArnaud 19961) . We model the SED with 
either a single or broken power-law intrinsic spectrum - for 
the latter we assume ySx - Po + 0.5 (Sari, Piran & Narayan 
1998). Extinction is modelled with the three commonly assumed 
extincti on curves of Milky Way, Large and Small Magellanic 
Cloud (iPei 1992h . Only the spectrum red-ward of Lycr line has 
been used in analysis. The broadband SED is best described by 
a broken power-law and an SMC-type extinction with Ay = 
0.11 + 0.02 mag and Po - 0.33 + 0.02. Requiring Ap - 0.5 
it has been possible to better constrain the X-ray column density 
to Nh - l-35;^[ Qg X 10^^ cm“^. This value is consistent with the 
direct estimate from X-ray data only (see Section 3.2). 

As a consistency check we also built the spectral energy dis¬ 
tribution of the optical afterglow using all the available photo¬ 
metric observations reported in Table |2] and Swift X-ray data. 
The extinction estimated with photometric information (Ay ~ 
0.6 mag) is much higher than the one obtained from the more 
accurate spectral analysis reported above. This is due to the 
fact that photometric measurements are too sparse to give re¬ 
liable results. The rest-frame extinction (Ay ~ 0.1 mag) of 
GRB 140515A is consistent with the Ay distribu tion found for 
the complete BAT6 sample (ICovino et al. 20131) and suggests 
that this is a typical value also for high-z events. 


4. Discussion 

4.1. Late time flare emission / refreshed shock 

The most straightforward explanation for the observed X-ray 
peak at t ~ 3000 s (Eig. |2]i is the onset of the afterglow emission. 
This interpretation is supported by the rise and decay indices of 
the X-ray light curve, that are consistent with the expectations 
for the emission of the forward shock interacting with an ho¬ 
mogeneous medium. An alternative possibility is that the X-ray 
peak corresponds to a late time flaring activity, or to variability of 
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^rest [10'® Hz] 


Fig. 5. Spectral energy distribution obtained with X-shooter 
(red) and XRT data (black). We included also the radio detection 
(green) of Laskar et al. (2014). Dashed lines mark the position 
of the injection frequency (v,- = 2 x 10*^ Hz) and the cooling 
frequency (v^ ~ 2 x lO'® Hz) expected for a pure synchrotron 
model. 


the GRB afterglow interacting with the ambient medium. A pow¬ 
erful tool to investigate the nature of this variability is the com¬ 
parison of the flux increase as a function of the temporal variabil¬ 
ity of the peak with the regions of allowance fo r bumps in the af - 
terglow on the basis of kinematic arguments (lloka et al. 2005h . 
In Fig. [6] we portrayed the s ample of early time {tp^ < 1 ks) 
flares (IChincarini et al. 2010l) and late time (tpk > 1 ks) flares 
dBernardini et al. 2011). 

If we interpret the broad X-ray bump of GRB 140515A as 
a single long-lasting flaring episode it would occupy a different 
region with respect to the observed X-ray flares since it is charac¬ 
terised by a very long duration (Af/fp*. ~ 50 » 1) and large flux 
variation (A/// ~ 10^). It would be t herefore consistent with 
being produced by refreshed shock s dRees & Meszaros 19981: 
iKurnar & Piran 2000bl : ISari & Meszaros 2000li 7 or by an intrin- 
sic angular structure on the emitting surface (a “patchy shell”) 
dMeszaros et al. 19^ iKumar & Piran 2000^ . or by shock re¬ 
flection generated by the interaction of the reverse shock 
with dense shells fo rmed at an earlier stage of the explosion 
dHascoet et al. 201"^ . An increase of the external medium den¬ 
sity would require a sharp and large jump in a uniform density 
profile to produce the observed increase in the observed light 
curves, which seems unlikely. 

A single X-ray broad peak is well outside the region of va¬ 
lidity for the internal shock model. However, there is still the 
possibility that the broad single peak that we are observing 
is the result of the superposition of multiple peaks, each with 
At It pic « 1. In the inset of Fig. |2] we sketched a possible tem¬ 
poral behaviour for GRB 140515A, where three flaring events 
(with a typical profile as described in iNorris et al. 20()5 ). super¬ 
posed to the underlying temporal decay, could be responsible 
of the shape of the broad bump observed. If we consider this 
scenario the observed behaviour becomes consistent with the in¬ 
ternal shocks scenario (Fig. This situation resembles the case 
of GRB 050904, a GRB at very similar redshift (z = 6.29) that 
shows a late time variability in the X-rays and a sudden drop of 
the observed emission afterwards. However, GRB 140515A is 



it/tp, 


Fig. 6. Kinematically allowed regions for afterglow variabil¬ 
ity in the A/// vs. Af/fpk plane. Coloured lines with ar¬ 
rows represent the allowed regions for density fluctuations 
on-axis (blue), density fluctuations off-axis (red), multiple 
density fluctuations off-axis (green), refreshe d shocks (pink) 
and patchy shell (black), respectively (se e lloka et al. 2005 : 
IChincarini et al. 2010l:lBernardini et al. 201 ll for details). In this 
plot we show early time {tpk < 1 ks, grey points) and late time 
{tpk > 1 ks, magenta squares) flares. The red triangles are the 
three flaring episodes in GRB 140515A. The error bars account 
for the uncertainty on the behaviour of the underlying contin¬ 
uum: the lower bar corresponds to a flat power-law decay after 
10^ s, the high bar to a flat decay normalised to the last datapoint, 
and the central value to a power-law decay consistent with the 
slope of the optical light curve after 10^ s. 


fainter than GRB 050904 and its variability has not been fully 
captured by the XRT. 

4.2. Standard afterglow interpretation 

In order to explain the observed light curves, we consider a 
semi-analytic model that describes the dynamical evolution 
of the fireball when interacting with the external circumburst 
medium and the respe ctive radiative emi ssion in the standard for¬ 
ward shock scenario (iNava et al. 2013h . The radiative descrip¬ 
tion is based on the model illustrated by Nappo et al. (2014), 
that allows us to compute the synchrotron spectrum as a func¬ 
tion of time, normalised to the bolometric luminosity obtained 
by the dynamical model. We assume that the electrons are in¬ 
jected with a power-law energetic distribution with index p and 
can cool for synchrotron and synchrotron self-Comtpon (SSC) 
radiation. The model allows to obtain at each step: 

(i) the synchrotron break frequencies (i.e. the self-absorption 
frequency Va, the injection frequency Vi and the cooling fre¬ 
quency Vc) 

(ii) the fraction of dissipated energy that is emitted in radiation 
erad that is used to determine the bolometric luminosity. 

(iii) the comptonization parameter Y 

(iv) the synchrotron spectrum Ly ^yn 

Since the spectrum of the radiation is estimated at each time, 
the light curve at a specific frequency Fy(tobs) can be derived. 
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Time [s] 

Fig. 7. Two multi-wavelength interpretations of GRB 140515A 
afterglow light curve in X-ray band (black dots, blue lines), 
SDSS-z' band (yellow circles), J band (orange circles), H band 
(red circles, red lines) and radio band (green star, green lines). 
Dashed lines represent the solution in the “flat electron spec¬ 
trum” scenario {p — 1.67). Solid lines describe the solution for 
the “multi-component” scenario {p - 2 . 1 ). 


Table 4. Parameters of the proposed light curves scenarios. We 
estimated Pjet from Equation 1 in Sari, Piran & Halpern (1999) 
and Ghirlanda, Ghisellini & Lazzati (2004), using Eiso = 6 x 
10 ^^ erg. 




flat spectrum 

multi-component 

To 


195 

125 

7 


0.02 

0.04 

n 

[cm 

0.05 

0.5 



0.15 

0.12 

fB 


2.3 X lO-'* 

1.5 X lO-'^ 

P 


1.67 

2.1 

^jet 

n 

>1.5 

>2.9 

Ey 

[erg] 

> 2.1 X 10-*® 

> 7.9 X 10'*“’ 


The parameters that can be varied in order to reproduce the 
observed light curves are: the initial bulk Lorentz factor (To), the 
isotropic prompt emitted energy Eiso, the prompt radiation efh- 
ciency ( 77 ), the circumburst density (n, in case of homogeneous 
medium), the injected electron spectral index (p), the fraction 
of dissipated energy distributed to the leptons (ee) and to the 
magnetic held (ee)- The model assumes an isotropic ejecta, so 
it can not reproduce geometrical features as jet break and side- 
expansion effect. 

From the modelling, two different interpretations of the 
multi-wavelength light curve of GRB 140515A are possible: a 
pure synchrotron emission by electrons with a very flat spectrum 
(with p < 2) or a multi-component model (with p > 2). We con¬ 
sidered only the radio, J, H and the X-ray band emissions, since, 
as shown in section 3.3, the optical SDSS-z' emission is strongly 
affected by absorption and the resulting flux is very likely under¬ 
estimated. 


4.2.1. A very flat injected spectrum 

We model the observations in the J, H, radio, and X-ray bands as 
synchrotron radiation produced in the forward shock. In order to 
obtain a successful modeling we need to assume that the electron 
injection spectrum is a power-law with a very flat index (p = 
1.67) that extends up to some maximum Lorentz factor jmax- 
Even in this case, the model cannot And a solution for the early 
time (very steep) X-ray emission. As suggested by e.g. Ghisellini 
et al. (2009) we assume that the X-ray light curve is composed 
by a late prompt component decaying with time as a power- 
law of slope > -3 (dominant for fob.s ^ 600 s) and by a second 
component interpreted as the actual X-ray afterglow emission 

(fobs ^ 600 s). 

A good description of the detections in J, H, radio, and X- 
ray bands (Fig.|2l dashed lines) is obtained using the parameters 
reported in Table |4] The bump observed in the X-ray light curve 
is due to the onset of the afterglow, corresponding to Eq = 195. 

The last detection in the X-ray band shows a sudden drop of 
the X-ray flux. The temporal index after the break is very steep 
(oblate = 3.9 + 0 . 6 , see table [TJ, and is not consistent with the¬ 
oretical predictions from jetted outflows. As an alternative we 
then suggest that this steepening is due to the passage of the 
maximum synchrotron frequency (associated to the maximum 
Lorentz factor of the electrons) in the X-ray band. To model 
the maximum frequency we assume that the ratio between the 
maximum and minimum Lorentz factors of the electrons is con¬ 
stant in time, and is equal to ymaxljmin — 2800. This interpre¬ 
tation allows us to fully describe the observations with no need 
for a jet break. Since the light curve does not show any break 
until 10 ^ s, this time represents a lower limit on the jet break 
time. The corresponding lower limit on the jet opening angle 
and on the collimation corrected energy Ey are reported in Table 
|4] T his lower limit on is consisted with the Ept - Ey correla¬ 
tion (iGhirlanda et al. 2004t) . Consistency with this correlation, in 
fact, requires a jet break at times larger than 3.5 x 10^ s, strength¬ 
ening the hypothesis that the steep break in the late X-ray light 
curve is not due to the jetted geometry of the outflow. 

4.2.2. Multi-component model 

An alternative interpretation could be obtained with a steeper 
electron injected spectrum. However, in this case it is not pos¬ 
sible to obtain a pure synchrotron solution that can explain cor¬ 
rectly the J, H, the X-ray emission and the radio detection. In 
fact, if we describe simultaneously the J, H and radio emissions, 
we underestimate the X-ray light-curve. As shown also in previ¬ 
ous sections, the peculiar shape of the X-ray light curve, charac¬ 
terised by an important time variability up to fobs ~ 10 “* s, sug¬ 
gests that the X-ray emission at those times could be probably 
caused by the composition of the standard afterglow emission 
in a forward shock scenario and some additional emission (for 
instance flares, or a long lasting prompt emission). 

Therefore, the time of the X-ray peak must be similar to the 
deceleration time computed by the model. Moreover, the pre¬ 
dicted X-ray emission cannot overcome the observed flux, since 
it must be a composition of the afterglow and additional emis¬ 
sions. In this scenario, we expect that the rise of the X-ray flux 
(at fobs ~ 10^ s) corresponds to the rise of the X-ray afterglow. 
Adding few constrains we obtain a compatible multi-wavelength 
prediction of the afterglow light curve (see solid lines in Fig.|7]l. 
The set of used parameters for this scenario is reported Table 0] 

In this scenario the X-ray afterglow prediction is below the 
observed data at early times, while the last observed detection 
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becomes compatible with the expected afterglow emission. This 
is consistent with the hypothesis of a multi-component X-ray 
emission, and we can assume that only for fobs > 2 x 10^ s we 
are observing a pure X-ray afterglow emission, not contaminated 
by long-lasting central engine activity. 

At the time of the Chandra observation, the X-ray afterglow 
flux predicted by this modelling is marginally consistent with the 
Chandra upper limit. A jet break around this time would make 
the energetic consistent with the prediction of the Epk - Ey cor¬ 
relation. However, since a break is not strictly required by our 
modelling, we conservatively set a lower limit on the jet break at 
fjet > 2 X 10^ s. The corresponding lower limits on the jet open¬ 
ing angle and on the collimation corrected energy are listed in 
table |4] 


4.3. Pop III or enriched Pop II progenitor 


GRB 140515A shows evidence of long lasting central-engine 
activity up to ~ 10“* s after the burst event. Its red- 

shift (z > 6) could suggest a Pop III star progenitor. 

These type of massive stars (M > IOOMq), that formed 
in the early universe at low metallicity (Z < 10“^), have 
been also proposed as pr ogenitor of the so- called ultra-long 
GRBs, i.e. GRB 111209A dOendre et al. 2013h. GRB 12 1027A 
(iHou et al. 2014 . and GRB 130925A dEvans et al. 20l4 . 


In this scenario, the long duration is the results of the time 
needed for the accretion and collapse mechanisms. In the hy¬ 
pothesis of such a GRB progenitor one should expect to detect 
a very low density environment with a density prohle domi¬ 
nated by the IGM. Another expectation for such massive col¬ 
lapsing stars is a long-lasting blackbody emission component in 
their spectra, wi th a typical aver age rest-frame temperature of 
IcTbb ~ 0.5 keV (iPiro et al. 20141) . This thermal emission would 
be in principle detectable by BAT and/or XRT if the redshift of 
the event is low. 


In the case of GRB 140515A observations do support the 
idea of a low density environment with negligible contribution 
from the host galaxy, but there are no hints for a particularly low 
value of the metallicity (see Section 3.4.3). Moreover, being at 
such a high-z we do not expect to detect the blackbody compo¬ 
nent with Swift instruments. Indeed we tested this possibility (see 
Section 3.1) but we did not hnd any improvement of the ht with 
the inclusion of a blackbody component in the prompt emission 
spectrum. Therefore, the hypothesis that GRB 140515A origi¬ 
nated from a Pop III star (or even from a Pop II star with envi¬ 
ronment enriched by Pop III stars) is unlikely. 


4.4. Reionization and escape fraction of ionizing radiation 

The distribution of intrinsic column densities of GRB hosts can 
be used to constrain th e average escape fraction of ionizing ra¬ 
diation from the hosts (IChen et al. 2007h . under the assumption 
that GRB sightlines, taken as an ensemble, sample random lines- 
of-sight from star forming regions in GRB hosts. At intermedi¬ 
ate redshifts (z > 2) the sample of GRB hosts from Chen et al. 
(2007) indicates that only in about 5% of all cases one expects a 
GRB sightline with log(AHi) < 18.5. With GRB 140515A be¬ 
ing only 1 out of 7 GRBS with z > 6 (and only 1 out of 3 
with measured HI column densities), it appears that high red- 
shift GRB hosts may have, on average, lower HI column densi¬ 
ties and, hence, higher escape fractions than their lower redshift 
counterparts. 


Table 5. Absorption properties of the GRBs with z > 5 (for 
the events marked with * the redshift was estimated photometri¬ 
cally). References: 1) Evans et al. 2010; 2) Perley et al. 2010; 3) 
Jakobsson et al. 2006; 4) Covino et al. 2013; 5) Salvaterra 2015; 
6) Hartoog et al. 2014; 7) Totani et al. 2006; 8) This work. 


GRB 


log(NHi) 

log(NH,x) 

Ay 

Ref. 



[cm--] 

[ 10"' cm-"] 

[mag] 


060522 

5.11 

- 

< 160 

- 

1 

071025 

<5.2* 

- 

49 ± 19 

<0.54 

1,2 

140304A 

5.283 

- 

< 120 

- 

1 

050814 

5.3 

- 

< 16.8 

<0.9 

3, 1 

131227A 

5.3 

- 

^20+220 

- 

1 

060927 

5.467 

- 

< 36 

< 0.17 

4, 1 

130606A 

5.913 

19.93 

< 30 

< 0.2 

5,6 

120521C 

6.0* 

- 

< 60 

<0.3 

5 

050904 

6.295 

21.6 

63*3^ 

0.15 ±0.07 

5,7 

140515A 

6.327 

< 18.5 

13 5+12.2 

0.11 ± 0.02 

8 

080913 

6.695 

19.84 

951? 

0.12 ±0.03 

5 

090423 

8.26 


102± 

<0.1 

5 

120923A 

8.5* 

- 

<720 

- 

5 

090429B 

9.4* 

- 

140 ± 10 

0.10 ±0.02 

5 



Fig. 8. Ay, Nh, and Nn/Ay ratio as a function of redshift. 
Black points are from Covino et al. (2013) for events with 
z < 4, while the remaining events (blue circles, purple stars) 
are listed in Table |5] GRB 140515A is marked with a red star. 
The solid/dashed gray lines in the middle panel represent the 
effect of the intervening material alo ng the line of sight (see 
ICamnana et al. 20151 ISalvaterra 2015h . 


More quantitatively, the Kolmogorov-Smirnov test for the 
two distributions of HI column densities - first from Chen et al. 
(2007) and the second of four z > 5.9 GRBs with measured Ahi 
values - shows that the two distributions are consistent with only 
9% probability. That probability raises to 30% if GRB 140515A 
is excluded. The importance of constraining the escape fractions 
in reionization sources is obvious, so a larger sample of z > 6 
GRBs with measured HI column densities would be highly de¬ 
sirable. 
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Such a sample would also serve as a direct test of reion¬ 
ization at z > 6, where constraints from high redshift quasars 
become scarce. A significant advantage of GRBs over quasars 
is in their low or negligible bias. While bright quasars, likely, 
do reside in the most massive, highly biased dark matter halos, 
GRBs hosts at high-z seem to sample the general galaxy popu¬ 
lation. Hence, constraints for the neutral hydrogen fraction ob¬ 
tained from the analysis of the IGM damping wing profile in the 
absorption spectra of GRB hosts can be expected to be more re¬ 
liable than the analogous constraints from the quasar proximity 
zones. 

In addition, constraints on the mean neutral fraction 
from observations of QSO proximity zones are, typically, 
lower limits (neutr al fraction can be larger if a quasar 
lifetime is longerl dBolton et al. 201 ll : [Robertson et al. 20131: 
[Robertson et al. 2015h . while constraints from GRBs are upper 
limits. Hence, the two observational probes are highly comple¬ 
mentary to each other (this is demonstrated by red and orange 
diamonds in Fig. Ell. 


4.5. High-z GRBs absorption properties 

In Table |5] we report the absorption properties - neutral hydro¬ 
gen column density in the host galaxy (Nhi), X-ray equivalent 
column density (Nh.x), and optical dust extinction (Ay) - for 
all known GRBs with redshift z > 5. As show in Table |5] host 
galaxies of high-z gamma-ray bursts have small but not zero ex¬ 
tinction, and the value of their optical extinction remains in fact 
rather constant (0.1 < Ay < 0.2 mag) over a broad range of red- 
shifts. 

It should be noted that no large values for Ay have been ever 
detected for GRBs at redshift > 5, and this is probably due to an 
observational bias, since it would be difficult to carry out optical 
follow-ups. Nevertheless, due to the general low metals abun¬ 
dances of the young galaxies at such a high-z, it is also plausible 
that large extinction are intrinsically less probable than at lower 
redshifts. 

While the Ay does not seem to evolve with redshift, there 
are no detected events with low Nh,x at z > 6. This effect can 
be naturally explained by the increase of absorption of the in¬ 
tervening systems along the line-of-sigh t ( Camnana et al. 201^ 
ICovino et al. 2013llCampana et al. 20151;ISalvaterra 2015 ). This 


mimics the evolution of the Ne/Ay ratio with redshift observed 
in Fig. [8] (bottom panel) up to z ~ 10. 


5. Conclusion 

We presented the multi-band spectroscopic and temporal analy¬ 
sis of the high-z GRB 140515A. The overall observed temporal 
properties of this burst, including the broad X-ray bump detected 
at late times, could be explained in the context of a standard af¬ 
terglow model, although this requires an unusually flat index of 
the electron energy spectrum (p - 1 .67). Another possible inter¬ 
pretation is to assume that an additional component (e.g. related 
to long-lasting central engine activity) is dominating the X-ray 
emission. In the latter case, the broad band observations can be 
explained using a more typical value of the spectral index for 
the injected electron spectrum (p = 2.1). Our modelling in this 
case shows that the central engine activity should cease at late 
times (~ 2 X 10^ s), when the X-ray afterglow starts to dom¬ 
inate the emission. In both scenarios the cooling frequency is 
expected to be between the optical and the X-ray energy bands 
(Vc ~ 2 X lO'® Hz) and the average rest-frame circum-burst ex¬ 
tinction (Ay ~ 0.1) resulted to be typical of high-z bursts. 


Our detailed spectral analysis provided a best estimate of 
the neutral hydrogen fraction of the IGM towards the burst of 
xh! < 0.002 and a conservative upper limit of the HI abun¬ 
dance in the GRB host galaxy of Nhi ^ 10'*'^ cm“^. These val¬ 
ues are slightly different from the ones estimated by Chornock 
et al. (2014b). In addition, the spectral absorption lines ob¬ 
served in our sp ectra are the weakest lines e ver observed in 
GRB afterglows (ide Ugarte Postigo et al. 2013) . suggesting that 
GRB 140515A happened in a very low density environment. 
However, our upper limits on the gas-phase abundances, cou¬ 
pled with the fact that we cannot establish the exact metal-to- 
dust ratio, do not allow us to distinguish between metallicity in 
the range of 10 < [Z/i/] <0.1. This makes the possible Pop 

III star origin for GRB 140515A uncertain and doubtful. 

For all high-z GRBs the contribution of the host galaxy was 
not negligible (TableQ. GRB 140515A is the first case when this 
does not happen, allowing us to give the the best observational 
constrains on a theoretical model at z > 6. 
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